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ABSTRACT 

Clusters of galaxies are mainly formed by merging of smaller structures, according 
to the standard cosmological scenario. If the mass of a substructure is > 10% of that 
of a galaxy cluster, the temperature distribution of the intracluster medium (ICM) in 
a merging cluster becomes inhomogencous. Various methods have been used to derive 
the two-dimensional projected temperature distribution of the ICM. However, meth- 
ods for studying temperature distribution along the line-of-sight through the cluster 
were absent. In this paper, we present the first measurement of the temperature stan- 
dard deviation along the line-of-sight, using as a reference case the multifrequency 
SZ measurements of the Bullet Cluster. We find that the value of the temperature 
standard deviation is high and equals to (10. 6± 3.8) keV in the Bullet Cluster. This 
result shows that the temperature distribution in the Bullet Cluster is strongly in- 
homogeneous along the line-of-sight and provides a new method for studying galaxy 
clusters in depth. 
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1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
structures with size ~1 3 Mpc (« (3 -5- 10)xl0 24 cm) 
containing a hot diffuse plasma (the intracluster medium or 
ICM) which eventually sets in equilibrium in the potential 
wells of the cluster. The plasma is hence heated to high tem- 
peratures during the cluster gravitational collapse and the 
subsequent virialization by the gravitational energy released 
by the cluster formation from the assembly of smaller struc- 
tures. The temperature of the ICM in its equilibrium stage 
can be as high as 10 -j- 15 keV in the most massive clusters 
and ensures that the elements present in the ICM are highly 
ionized. Abundant light elements (such as hydrogen, helium, 
and oxygen) in the ICM have all the electrons removed from 
their nuclei. The free hot electrons in the ICM cause X-ray 
bremsstrahlung emission due to the interaction with protons 
(see, e.g., Rybicki & Lightman 1979) and also cause a dis- 
tortion of the CMB radiation towards galaxy clusters due 
to the Inverse Compton (IC) scattering off CMB photons 
(the Sunyaev-Zel'dovich effect, hereafter the SZ effect; see, 
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e.g., Sunyaev-Zel'dovich 1980). The typical number density 
of free electrons in the ICM is 10~ 3 -f- 10~ 2 cm' 3 . 



Mergers of sub-clusters occurring during cluster evolu- 
tion cause an additional significant and non-uniform heat- 
ing of the ICM due to the conversion of kinetic energy of 
gravitationally accelerated sub-clusters into thermal energy 
released into the ICM by shock waves. Shocks have been 
found, in fact, by the Chandra X-ray observatory in var- 
ious systems, from galaxy groups with temperature of ~1 
keV to the most massive galaxy clusters with temperature 
of «15 keV (for a review, see Markevitch & Vikhlinin 2007). 
The cluster 1ES0657-558 (known as the Bullet Cluster) is a 
"textbook" example of a hot, merging cluster which has a 
bullet-like substructure in the ambient ICM, a strong shock 
wave driven by the bullet substructure, and very high and 
non uniform values of temperature (Markevitch et al. 2002). 
Thanks to deep Chandra observations, the Bullet Cluster 
is one of the best-studied galaxy clusters in X-rays. The 
analysis of the X-ray data from Chandra (see, e.g., Million 
& Allen 2008) shows that a single temperature model (i.e., 
a single temperature value observed along the line-of-sight 
through the cluster) does not properly describe the observed 
its X-ray spectrum. The temperature distribution along the 
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line-of-sight in this cluster is then presumably strongly in- 
homogeneous. 

Gas temperature in galaxy clusters can be derived by 
using different techniques, such as analysing of the X-ray 
spectra or the spectral form of the CMB distortion owing to 
the SZ effect. The derived X-ray spectroscopic temperature 
for a multitemperature electron population is weighted with 
a factor that is proportional to the squared number density 
and to the gas temperature to the power fs-3/4 (Mazzotta 
et al. 2004). Therefore, due to the non-linear combination 
of density and temperature there is no simple way to in- 
terpret the X-ray spectroscopic temperature as the actual 
gas temperature and to obtain direct information about the 
temperature distribution along the line-of-sight. Additional 
assumptions about the gas density and temperature pro- 
files along the line-of-sight, which are a priori unknown in 
merging galaxy clusters, are required for the X-ray analy- 
sis (see, e.g., Samsing et al. 2012). The temperature derived 
through the SZ technique is the mass-weighted temperature 
(Hansen 2004; Prokhorov et al. 2010, see also Colafrancesco 
& Marchegiani 2010 for a discussion on the spectroscopic 
SZ temperature derivation) and, therefore, it has a more ex- 
plicit interpretation. Below, we will use this fact to obtain 
the information about the temperature distribution along 
the line-of-sight. 

The SZ effect is an important tool for studying of galaxy 
clusters, since both the amplitude and spectral form of the 
CMB distortion (caused by the IC scattering of CMB pho- 
tons off free electrons) depend on the ICM physical pa- 
rameters (e.g., the plasma pressure and temperature). Rel- 
ativistic corrections to the SZ effect are significant for high- 
temperature plasmas in massive clusters of galaxies (see 
Rephaeli 1995). The contribution from the relativistic cor- 
rections becomes stronger at high frequencies. The first de- 
tection of the SZ effect brightness increment from the Bul- 
let Cluster at frequencies above 600 GHz has been obtained 
with the Herschel-SPIRE instrument (Zemcov et al. 2010). 
These high-frequency SZ observations have been used by Co- 
lafrancesco et al. (2011) to demonstrate that plasma models 
with several electron components provide a better fit than 
that with a single thermal component. This means that a 
substantial temperature variance (or standard deviation) 
a is expected along the line of sight through this cluster. 
The method of measurements of the temperature standard 
deviation, a, along the line-of-sight has been proposed by 
Prokhorov et al. (2011). The use of the temperature stan- 
dard deviation, which is a measure of the plasma inhomo- 
geneity along the line of sight, provides us hence with a 
model-independent approach to characterize the amount of 
the temperature inhomogeneity in the ICM. 

In this paper, we present the first measurement of the 
temperature standard deviation a along the line-of-sight in 
galaxy clusters. We use the method of Prokhorov et al. 
(2011) and SZ effect intensity data for the Bullet Cluster 
measured at four frequencies (150, 275, 600, and 857 GHz) 
to derive the temperature standard deviation. Using this 
method, we demonstrate that the temperature distribution 
along the line of sight through this cluster is very inhomo- 
geneous and we quantify such an inhomogeneuty by using 
the value of the standard temperature deviation. 



2 THE TEMPERATURE STANDARD 

DEVIATION OF THE BULLET CLUSTER 

In this section, we calculate the temperature standard devi- 
ation for the Bullet Cluster. This high temperature, merging 
galaxy cluster is a good target for measuring the standard 
temperature deviation because the contributions of high or- 
der corrections to the SZ effect are strong and because the 
gas temperature distribution is expected to be inhomoge- 
neous along the line of sight owing to the cluster merger 
effects. In the following, we first describe the method to cal- 
culate the temperature standard deviation along the line of 
sight that takes into account SZ measurements at different 
frequencies and then apply this method to the Bullet Clus- 
ter. 

The relativistically correct SZ effect can be described 
in the formalism based on the extension of the Kompa- 
neets equation which allows relativisitic effects to be in- 
cluded (e.g., Challinor & Lasenby 1998; Itoh et al. 1998). 
Analytical forms for the spectral changes due to the SZ ef- 
fect, that are correct to second order in the expansion pa- 
rameter O = fc b T c /(m c 2 ), can be written as 



AI(x) 
Io 



' (Jt 



J n c ex(j (i)+ex 9 i(i)+e 2 xj 2 (i))iii, (1) 



where Io = 2(fc b T cm b) /(he) , x = hv/k h T crnh , T c is the 
electron temperature, n gas is the number density of the ICM, 
or is the Thomson cross-section, m e the electron mass, c the 
speed of light, fcb the Boltzmann constant, and h the Planck 
constant. The spectral functions of go(x), gi(x), and g2(x) 
are taken from Eqs. (28) and (33) of Challinor & Lasenby 
(1998) and are given by 
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(4) 



where C(x) = x x coth(a;/2). 

The spectral function, go(x), corresponds to the spec- 
tral function derived from the Kompaneets approxima- 
tion (Zel'dovich & Sunyaev 1969). Both the first-order and 
second-order relativistic effects make significant contribu- 
tions to the spectral distortion for k h T c > lOkeV. We veri- 
fied that, for a plasma with the temperature of khT c = 13.9 
keV and optical depth r = 1.3xl0~ 2 , as used by Zemcov et 
al. (2010), the second-order corrections in r to the spectral 
shape (see Colafrancesco et al. 2003, for technical details) do 
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not produce appreciable variations in the SZ spectrum and, 
therefore, their contribution will be neglected in our paper. 

The CMB intensity change produced by the SZ effect 
from a galaxy cluster with inhomogeneous density and tem- 
perature distributions in the formalism based on an exten- 
sion of the Kompaneets equation, that is correct to second 
order in the expansion parameter O, is given by 

A/^) f(k b T e ) ((k b T c ) 2 } 

— « r ± {- x g {x) + 11 x gi(x)+ 

Io \ m c c z mid 1 

\^x 52 (x)j, (5) 

where (khT c ) = / n c khT c dl/ j n c dl is the tempera- 
ture averaged along the line-of-sight and {(k^T c ) 2 } = 
J n e (kbT e ) 2 dl j J n e dl is the squared temperature averaged 
along the line-of-sight. Note that the optical depth of the 
electron plasma, r, can be derived from X-ray observations. 

To study the gas inhomogeneity along the line of sight, 
we calculate the temperature standard deviation, a, us- 
ing the method based on SZ intensity measurements (see 
Prokhorov et al. 2011). 

The temperature standard deviation defined as 

o = V((k b T c ) 2 ) - ((k h T c )) 2 (6) 

can be derived from the first two terms in Eq. ((5J by using 
the SZ measurements at three or more frequencies. The value 
of this quantity for a homogeneous gas temperature distri- 
bution equals zero by definition. A finite, non-zero value of 
a is therefore a measure of the inhomogeneity of the gas 
temperature along the line of sight. 

We derive the value of the temperature standard devi- 
ation by using the available multifrequency SZ observations 
of the Bullet Cluster: these have been obtained by ACBAR 
at frequencies of 150 GHz and 275 GHz (Gomez et al. 2004) 
and by Herschel-SPIRE at frequencies of 600 GHz and 857 
GHz (Zemcov et al. 2011). This set of SZ observations have 
been used by Colafrancesco et al. (2011) in their studies of 
plasma models with different thermal and non-thermal com- 
ponents. The central SZ intensity values, Al /Io, measured 
by ACBAR and Herschel-SPIRE are shown with errorbars 
in Fig. [T]and listed in Table[T](in units of MJy sr _1 ). Note 
that the central SZ intensity values have been derived by 
using the radial profiles with an isothermal /3 model (for a 
review, see Birkinshaw 1999) taking the core radius and ft 
parameter from Tucker et al. (1998) for the ACBAR data 
(Gomez et al. 2004) , and the density profile given by Halver- 
son et al. (2009) for the analysis of the Hershel-SPIRE data 
(Zemcov et al. 2010). We checked the consistence between 
the central SZ intensity values derived from these models 
and found that the difference in the central values is ~ 10%. 
However, the gas temperature distribution should be inho- 
mogeneous along the line of sight and isothermal /3 models 
are just an approximation used to fit the data. With that 
caveat in mind, we will take the central SZ intensities from 
the papers by Gomez et al. (2004) and Zemcov et al. (2010) 
as it has been done in Colafrancesco et al. (2011). 

The SZ intensity at each frequency is given by the sum 
of zero-th, first, and second order contributions in <d and 
can be written as 

AI(x-) 

— pzAqx g (xi)+Ai x gi(xi) + A 2 x g 2 (xi) (7) 
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Figure 1. The SZ spectrum of the Bullet Cluster fitted with 
the zero-order term (dotted), with the zero plus first order terms 
(dashed), and with the zero, first, and second order terms (solid). 



Table 1. The central values of the SZ effect at the various fre- 
quencies as obtained by ACBAR and HERCSHEL-SPIRE 



Frequency (GHz) 


A 7 (MJy sr" 1 ) 


Error (MJy sr -1 ) 


150 


-0.325 


0.015 


275 


0.21 


0.077 


600 


0.268 


0.031 


857 


0.097 


0.019 



where the dimensionless frequencies x\, X2, x^, and £4 
equal to 2.64 (^=150 GHz), 4.84 (^=275 GHz), 10.57 
(^=600 GHz), and 15.09 (^=857 GHz), respectively. We 
use the function "MPFITEXPR'Q to find the best set of 
model parameters (Aq, Ai, and A2) which match the data 
performing Levenberg-Marquardt least-squares fit. We also 
apply the same method to find the best set of parameters 
in the models based on the Kompaneets approximation and 
on the expression including the zero and first order terms 
from Eq.©, respectively. We denote these models as Ml, 
M2, and M3, accordingly to the used number of terms pro- 
portional to O for calculating the SZ effect (i.e., the model 
Ml corresponds to the Kompaneets approximation). The 
best-fit sets of parameters (Aq, Ai, and A2) and 1-sigma 
errors on these parameters (SAo, SAi, and 8A2 ) derived 
from the SZ data are shown in Table 2 for the models of 
Ml, M2, and M3. The values of x 2 for these models used to 
fit the observed SZ spectrum of the Bullet Cluster are also 
shown in Table 2. Comparing the values of \ 2 from Table 2 
with those from the x 2 table, we conclude that the model 
M3 agrees with the observed SZ data much better than 
the models Ml and M2. The model M3 is consistent with 
the SZ data and the \ 2 value in this case is 0.3706 with 1 
degree of freedom (d.o.f.). The SZ spectra derived from the 
models Ml, M2, and M3 are shown in Fig. Q] by the dotted, 
dashed, and solid lines, respectively. Comparing the SZ 
spectra derived from the models with the observed SZ data, 
we find that the high-frequency observations provided by 

1 http : //cow. physics .wise . edu/~ craigm/idl/idl .html 
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Table 2. The values of parameters (Ao, A\, and Ai), errors on 
these parameters, and \ 2 for the models Ml, M2, and M3 used 
to fit the observed SZ spectrum. 



Parameter 


Ml 


M2 


M3 


A 


3.338965E-4 


3.66015E-4 


3.68631E-4 


Ai 




1.05000E-5 


1.46551E-5 


A 2 






5.21968E-7 


x 2 


30.39981 


8.42973 


0.37056 


d.o.f. 


3 


2 


1 


SA 


1.45174E-5 


1.56270E-5 


1.56461E-5 


SAi 




2.23969E-6 


2.6751E-6 


8A 2 






1.8488-7 



Herschel- SPIRE allow us to constain the contributions of 
the first-order and second-order relativistic SZ corrections 
to the thermal SZ effect. We calculate the covariance matrix 
for the parameters, Ao, Ai, and Ai, in the model M3 using 
the function of "MPFITEXPR" . The calculated covariance 
matrix is 

/2.44800S - 10 1.42381S - 11 1.67978S - 13\ 
co Vi j= 1.42381S - 11 7.15626S - 12 2.68991S - 13 
\1.67978£- 13 2.68991B- 13 3.38113.E - 14/ 

where i, j = 0, 1, and 2. The square root of the diag- 
onal elements of the covariant matrix gives the 1-sigma 
statistic errors on the parameters Ao, Ai, and Ai. 

Using the definitions of the parameters of Ao and Ai, 
we rewrite Eq.([5| as 



m c c 



At 



(8) 



We derive the value of the optical depth, r, using the 
X-ray observations of the Bullet Cluster by ROSAT those 
have been analyzed by Tucker et al. (1998). The surface 
brightness profile of the primary cluster was obtained by 
extracting the net counts in concentric annuli about the 
primary X-ray peak after excluding the western quadrant 
(which contains the secondary peak) and all point sources. 
Tucker et al. (1998) fitted the observed X-ray surface bright- 
ness profile to the standard hydrostatic-isothermal j3 model 
and obtained the values of a core radius, a central electron 
number density, and a beta parameter for the best fit. We 
re-scale the values of a core radius, a central electron num- 
ber density using the present-day Hubble constant value, 
Ho — 74 km s _1 Mpc -1 . Using the parameters of the beta 
model and Monte-Carlo simulations, we find that the value 
of the optical depth is r = 0.0138 ± 0.0016. 

To calculate the temperature standard deviation along 
the line of sight in the Bullet Cluster, we use the derived 
values of parameters (Ao, A\, and A 2 ), 1-sigma errors on 
these parameters, covariance matrix (cov), and optical depth 
(r = 0.0138 ± 0.0016). The 1-sigma error on the standard 
temperature deviation is given by 



5a 



\ i=0,l 



9it\ ... „ da da ( 9cr\ , 

dTj 5A ' + 2 8A -0A 1 COVl2+ {^ ] - 



(9) 

We consider the value of the optical depth as derived from 



X-ray observations, while the parameters (Ao, Ai, and A 2 ) 
are derived from SZ observations. Thus, the optical depth is 
an uncorrelated variable with other parameters in Eq.©. 

Using Eqs.© and ((9j), we find that the value of tem- 
perature standard deviation in the Bullet Cluster equals to 
a = 9.5 ± 2.6 keV. This value of standard temperature de- 
viation is high although somewhat lower than the average 
plasma temperature, k^T c « 14.5 keV, derived from the 
Chandra X-ray observations of the Bullet Cluster (see, e.g., 
Million & Allen 2009). 

Since the derived value of temperature standard devi- 
ation is high and the contribution of higher order terms in 
the expansion parameter O is significant at high frequencies 
for temperatures ^ 15 keV, we recalculate the values of 
the temperature standard deviation and 1-sigma error on 
this quantity considering the spectral changes due to the 
SZ effect correct to third order in O. We estimate the un- 
certainty of the technique comparing the SZ intensity value 
derived in the extended Kompaneets formalism with that is 
derived from the Wright formalism (1979) and take them 
into account as systematic errors in the analysis. In this 
analysis, we consider electron temperatures in the range 5 
keV - 30 keV that agree with that was found in the recent 
hydrodynamic simulations of the Bullet Cluster (Akahori & 
Yoshikawa 2011). We find that the values of the temperature 
standard deviation and 1-sigma error, a = 10.6 ± 3.8 keV, 
are consistent with those derived above. 

We also calculated the temperature standard deviation 
using the optical depth derived by using the radial density 
profiles given by Halverson et al. (2009) and Ota & Mitsuda 
(2004) . We obtained the values of the temperature standard 
deviation and of the associated errors that are within 10% 
and 30%, respectively, of those derived above. Therefore, 
the temperature standard deviation is high compared to the 
value of the gas temperature also assuming these models for 
the radial density distribution. 

Based on our results, we conclude that the plasma tem- 
perature distribution is very inhomogeneous along the line 
of sight in the Bullet Cluster. This result confirms the pres- 
ence of additional electron components previously revealed 
by Million & Allen (2009) by means of X-ray observations 
and by Colafrancesco et al. (2011) by means of SZ observa- 
tions. Note that the method of Prokhorov et al. (2011), that 
we use in this paper, can be also applied to other merging 
massive galaxy clusters as soon as high-frequency SZ obser- 
vations of other clusters become available. 



3 CONCLUSIONS 

We have derived the value of temperature standard devia- 
tion, a, along the line of sight in the Bullet Cluster using SZ 
effect observations at frequencies of 150, 275, 600, and 857 
GHz. The measured value, a = 9.5 ± 2.6 keV, is obtained 
from the analysis taking into account the three first terms 
in the expansion parameter O = ki,T c /mc 2 of the SZ effect 
(see Eq. [1]), and the analysis including the four first terms in 
O and taking into account the uncertainty of the technique 
gives a consistent value, a = 10.6 ± 3.8 keV. The derived 
value of a is of the order of the average plasma tempera- 
ture, kbT e = 14.5 keV, in this cluster. Such a high value of 
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a suggests hence that the gas temperature distribution is 
strongly inhomogeneous along the line of sight. 

The high plasma temperature inhomogeneity along the 
line of sight in the Bullet Cluster is likely caused by its past 
merger activity. The presence of the strong shock wave and 
"bullet" cold substructure in this galaxy cluster is strong 
evidence in favor of a recent major merger of subclusters. 
Our result shows that the value of temperature standard 
deviation along the line of sight can be quite high in mas- 
sive clusters undergoing strong merging and, therefore, that 
the plasma in such clusters should be strongly disturbed. 
This agrees with the fact that a single temperature model 
does not provide a good fit to X-ray observations (see Mil- 
lion & Allen 2009) and to SZ effect measurements (see Co- 
lafrancesco et al. 2011). By comparing the derived temper- 
ature standard deviation for the different models of a radial 
density distribution, we verified that our results are insen- 
sitive to the choice of the radial density profile. However, 
note that our results can be subject to additional possible 
systematic errors that come from: i) the uncertainties in 
the flux scales between the different frequency bands; ii) in- 
complete knowledge of the radial gas density distribution; 
iii) the subtraction of the bright sub-mm sources that can 
cause confusion in the ACBAR as well as the SPIRE data 
(see Johansson et al. 2010). The first two of these issues will 
be considered in a following paper (Prokhorov et al. 2012) 
where a joint analysis of the X-ray and SZ data for the Bullet 
Cluster will be presented. 

The technique we present here does not require any spe- 
cific assumptions about temperature profile. Therefore, the 
advantage of using the temperature standard deviation is 
that this quantity provides a model-independent informa- 
tion of the plasma temperature inhomogeneity along the line 
of sight. 

There are several astrophysical and cosmological con- 
sequences of having large values of temperature standard 
deviation, a, in galaxy clusters: 

i) measuring a in different clusters of galaxies will tell us the 
relevance of merging and shocks in the evolution of these 
systems, especially those undergoing a collision along the 
line-of-sight; 

ii) the specific value of a will allow us to constrain the pos- 
sible non- gravitational effects in clusters and their evolution 
with the mass and redshift of the system. In addition it will 
tell us how much of the energy density of the cluster is due to 
non- gravitational effects with respect to pure gravitational 
effects; 

iii) the value of a will tell us how confidently we can use 
galaxy clusters as cosmological probes (e.g., those are based 
on using of a mass-temperature relation) assuming simple 
models for the ICM, such a single temperature model. 

We will discuss these and other consequences in more details 
elsewhere. 

The first measurement of temperature standard devia- 
tion along the line-of-sight, that we have presented in this 
paper, demonstrates that existing microwave and mm. in- 
struments provide us with a unique opportunity to derive 
this quantity in massive, merging clusters. Future multi- 
frequency SZ effect observations including high frequency 
observations of other merging galaxy clusters will allow us 
to study the plasma temperature distribution along the line 
of sight for many clusters of galaxies and will provide us with 



a better understanding of gas dynamical processes occurring 
in the ICM of merging clusters. 

This is the first attempt in which we have been suc- 
cessful to measure a, and we plan to provide a systematic 
analysis of the temperature standard deviation in all clusters 
with sensitive enough multi-frequency SZE observations. We 
believe that this new tool will provide a promising method 
to analyze the complex structure of cluster atmospheres. 
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